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INTRODUCTION 


f The analysis of exposed surfaces on LDEF since its retrieval in 1 990 has revealed a 
I wide range of meteoroid and debris (M&D) impact features in the sub-micron to 
I millimeter size range, ranging from quasi-infinite target cratering in LDEF metallic 
structural members ( e.g . inter-costals, tray clamps, etc.) to non-marginal perforations in 
metallic experimental surfaces (e.g. thin foil detectors, etc.). Approximately 34,000 
impact features are estimated to exist on the exposed surfaces of LDEF. The vast 
majority of impact craters in metal substrates exhibit circular footprints, with -50% 
^retaining impactor residues in varying states of shock processing. 

Four papers (Berthoud et al., 1993; Hill etal., 1993; MacKay et at. , 1993; Paul, 
1 993) at the last (3rd) LDEF Post- Retrieval Symposium indicated that understanding of 
crater geometry and impactor residues as functions of impact velocity and angle can aid 
in the derivation of impactor trajectories and thus the directional distribution of M&D 
particles. 

The fundamental goals of this project were to duplicate and analyze meteoroid 
impact damage on spacecraft metallic materials with a view to quantifying the residue 
retention and oblique impact morphology characteristics. Using the hypervelocity impact 
test facility established at Auburn University a series of impact tests (normal and oblique 
incidence) were executed producing consistently high (11-12 km/s) peak impact 
velocities, the results of which were subsequently analyzed using Scanning Electron 
Microscope (SEM) and Energy Dispersive X-ray Spectroscopy (EDXS) facilities at 
Auburn University. 

2. OBJECTIVES 

The objectives of this program of research have been twofold: 

• to simulate and characterize meteoroid and debris hypervelocity impacts 

• to engineer and upgrade the Auburn University hypervelocity impact facility. 

The major issue for the LDEF meteoroid and debris investigators was that the impact 
angle is unknown and therefore only a coarse estimate of M&D flux directionality was 
possible. For impact velocities in excess of -7 km/s (the vast majority) there is very little 
residual evidence of impact angle effects, i.e. craters become progressively symmetrical 
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with no elongation along.the impact axis as the impact velocity increases. Figures 1 and 
2 show impact craters produced by oblique hypervelocity impacts by olivine particles 
(diameter = 20- 100 pm) onto copper targets. Most of the hypervelocity impact sites on 
the LDEF spacecraft occurred on quasi-infinite metallic surfaces and therefore this 
project has focused on such targets. Copper was selected as the target material to allow 
the separation of Mg, Fe elemental residue components from the Cu substrate and the 
identification of Al contamination from the hypervelocity impact facility drive process. 

One way to extract some further directional information is to analyze the impactor 
residues for any asymmetry in distribution. Figures 3 and 4 show the potential residue 
distributional effects associated with an oblique hypervelocity impact. The hypothesis is 
proposed that there will be some degree of quantifiable translation of residue material 
downrange for oblique impacts which should correlate with the impact velocity and 
impact angle. 

3. EXPERIMENTS 

3.1. Hypervelocity Impact Phenomena 

When a hypervelocity particle impacts a solid surface and intense flash of light is 
produced. Early researchers concluded that the impact-induced flash intensity increased 
with impact velocity, experimental data revealing a functional relationship between 
velocity, mass and flash intensity (Clark et al., 1959; MacCormack, 1963; Gehring & 
Warnica, 1963). 

The flash comprises two sub-components: (i) a short-duration transient spike 
attributable to a high pressure, high temperature, dense plasma, with resonance 
(collision) broadening of line structures, at the impact point, and (ii) a slow-rising, long- 
duration tail, exhibiting no line-broadening, due to the radiation of a neutral gas 

i 

expanding from the impact zone (Jean & Rollins, 1970). Observation of the impact flash 
in the experiments presented here allows the detection of particle impact location and 
the determination of impact velocity. 

The phenomenological model of hypervelocity impact was summarized in the 
mid-1 960’s (Gehring etal., 1965), based on the shock structures induced in the impactor 
and target. Essentially, the impact velocity (in combination with the impactor/target 
materials) determines the shock compression. The material of the impactor and target 
are compressed along the shock Hugoniot, increasing entropy, but are released 
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adiabatically, resulting in a residual temperature of between -10,000 - 50,000 K on final 
release which is sufficient to totally vaporize and partially ionize the impactor material 
plus a significant portion of the material excavated from the impact crater. Material that 
is not ejected from the impact crater will cool on release and will coat the interior of the 
crater, including the raised outer crater lips. 

3.2. Apparatus 

3.2.1 . Hypervelocity Launcher 

For the experiments reported here the hypervelocity impact facility at the Space Power 
Institute, Auburn University, Alabama was used. It utilizes the thermal expansion of a 
plasma produced by an electrically-exploded aluminum foil to accelerate micro-particles 
in the size range 10 p.m - 400 urn diameter. The system is described in detail elsewhere 
along with its typical operational sequence (Rose et at., 1991). A schematic of the 
system is shown in Figure 5. The plasma exit velocity at the muzzle has been measured 
at velocities in excess of 40 km/s using high-speed photographic techniques (IMACON 
790 image converter camera). The resultant velocity distribution of micro-particles is 
dependent on the combination of stored energy, particle size and density, and the 
number of particles under acceleration. 

The launcher is housed in a vacuum chamber pumped to a pressure of -10' 6 
torr. The micro-particles and residual launcher-drive gas pass along a 7 m long flight 
tube to a large (1 .45 m x 2 m) evacuated (~10‘ 6 torr) chamber where they are allowed to 
impact the target specimen or system under test. Electrical isolation is achieved since 
all gun electrical activity subsides before particle impact occurs in the target chamber. 

Furthermore, the hot plasma drive has relaxed and cooled and is captured by a 
series of baffles placed in the flight line. As a result, cleanliness is exceptional for a 
launcher system of this type. A final pre-impact baffle with an aperture correctly sized 
and positioned to control the impact site location may be inserted into the flight line. 

3.2.2. Impact Location and Velocity Diagnostics 

The primary diagnostic for the hypervelocity impact experiments is a Hadland-Photonics 
IMACON 790 Image Converter Camera. The IMACON is utilized in streak mode to 
detect the incident micro-particles velocities and locations by observing the luminous 
impact-induced plasma plume. The plasma plume is detected by the IMACON, looking 
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parallel to the target surface, via a slit that limits the field-of-view to no more than 100 pm 
at the target surface. 

Imaging is triggered to commence at the time of launch of the particles and 
terminates -0.5 ms after the arrival of the last hypervelocity particles. The arrival of 
each hypervelocity particle produces a light flash on the photographic image with a 
horizontal (time-axis) displacement proportional to the impact velocity. Since the flight 
distance is known to be 7.05±0.05 m the velocity may be calculated from the time-of- 
flight. 


Impact site location is derived from the vertical displacement of the flash on the 
streaked photographic image of the slit (see Figure 6). Each impact produces two flash 
images one displaced positively, the other displaced negatively, from the image 
horizontal center-line. Light emitting diodes (LEDs) are placed at known positions such 
that their continuous images provide boundaries on the photographic image 
corresponding to the impact zone limits. 

3.3. Experiments 

Each experimental target comprised nine (9) polished copper target plates, each 
measuring 1.5" X 1.5°, arrayed in a in a 3X3 configuration. One set of nine plates was 
exposed to the hypervelocity particle flux in the impact facility target chamber per 
experimental shot. On completion of the impact experiment the target was removed 
from the target chamber and the impact sites located optically, catalogued, and 
correlated with velocity data derived from the impact experiment streak record. 

The plates were then removed from the target fixture and inserted into an SEM 
for further analysis. Although each site was imaged in real-time, it was not necessary to 
produce a photographic record of each site. For the purposes of this project it was 
necessary to obtain quantitative EDX data for each site. Early experiments focused on 
only two locations within the craters, but as the analysis technique became more 
sophisticated five locations per crater were analyzed. Impact experiments were 
repeated to cover a variety of impact angles varying from 30° to 75° relative to the target 
surface normal vector. Figure 7 shows the analysis site locations for the final set of 
impact craters analyzed in this project. 

Olivine particles (75 pm nominal diameter), having the chemical composition 
(Mg:Fe) 2 Si 04 and a Mg:Fe:Si ratio of -19:1:10, were launched at copper targets. 
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Olivine, although higher in density than that usually accepted for meteoroids (p = 0.5-1. 5 
g.cnrr 3 ), was chosen as the best meteoroid simulant material that can be accelerated in 
the hypervelocity impact launcher system. Note that the mass density of olivine is 3.21 
g.cnrr 3 and that of copper is 8.92 g.cnr 3 , resulting in a mass density ratio of -0.36. This 
figure corresponds to that for low density (p - 1 g.cnr 3 ) meteoroids impacting aluminum 
(p -2.7 g.cnr 3 ) targets. It should be remembered that the density ratio scaling between 
different materials combinations is a first order approximation only. Peak impact 
pressure, post-shock mass velocity, and shock velocity all play roles that will lead to 
deviations from density ratio scaling at second order level. 

A thin Mylar® film (thickness 0.5 pm) was placed in front of the target plates, 
such that it was always normal to the particle flight path. The hypervelocity particles 
pass through the film before striking the plate, allowing the particle size to be found 
(Rose et al., 1992). The film also minimizes gun debris contamination of the target 
surface. 

4. TECHNICAL ANALYSIS 

4.1 . Meteoroid and Debris Environments 

The space particulate environment in low Earth orbit (LEO) comprises two components; 
the interplanetary dust particulate (I DP) component and the space debris particulate 
(SDP) component. Two standard models of the IDP component exist: the NASA SP- 
8013 model (Cour-Palais etal., 1969) and the ESA-Grun model (Grun etal., 1985). The 
accepted model for the SDP component is the NASA TM-1 00471 model (Kessler et al., 
1988). 


For the purposes of this project it was critical to know the impact velocity and 
impact angle of the SDP and IDP components at any specified spacecraft surface. 
Impact velocity, combined with impactor and target material type (density and porosity) 
determines the strength of the shock waves propagating both forward into the target and 
rearward into the impactor. Grun et al. (1985) determined that the average velocity for 
an IDP entering terrestrial space is 20 km/s, resulting in average impact velocities (for 
normal incidence) on a spacecraft RAM surface of between 27.8 km/s and 27.4 km/s for 
the altitude range 250 km - 850 km. The corresponding average impact velocity range 
(for normal incidence on the WAKE surface) is between 12.2 km/s and 12.8 km/s. 
Recently, McDonnell (1992) determined the mean normally-resolved impact velocity of 
IDPs incident on the LDEF spacecraft to be 21.5 km/s for the RAM (East) surface and 
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11.1 km/s for the WAKE (West) surface. A value of 15.7 km/s was calculated for the 
North-, South-, and Space-facing surfaces. 

For the space debris component, Kessler (1984) computed an impact velocity at 
500 km (for 30° and 60° inclination orbits) distribution by analyzing the collision cross- 
sections for all intersecting debris orbits. The worst-case environment is for polar orbits 
(inc - 90°). The average impact velocity of space debris for a spacecraft in a 28.5° 
inclination orbit is 9.8 km/s, increasing to 1 2.5 km/s for a polar orbit. At present, there 
are no data that verify these models in terms of absolute magnitudes, although it is 
reasonable to assume that the relative magnitudes are valid based on the orbital 
mechanics of the scenario. 

Recently, Kessler (1992) proposed the existence of small-scale (1-100 pm 
diameter) space debris in highly elliptic orbits (HEOs) to account for the 13% of impact 
sites on the WAKE surface of LDEF being of terrestrial origin (Bernhard et ai, 1993). 
Such SDPs are attributed to break-ups of upper stages in geo-synchronous transfer 
orbits (GTOs) (perigee altitude = 250 km, apogee altitude = 35900 km ). The debris 
produced then cascades down through lower orbits under the combined effects of solar 
radiation (at apogee) and atmospherics drag (at perigee), resulting in SDP impact 
velocities in excess of those determined for primarily circular orbits. 

A detailed computer model of the IDP flux component was implemented to allow 
the computation of impact velocity and impact angle distributions for both RAM and 
WAKE surfaces. Both the NASA and ESA IDP models assume isotropic flux 
distributions with a geocentric velocity distribution based on the Erickson meteor velocity 
distribution (Erickson, 1968) (see Figure 8). The model was implemented for a gravity- 
gradient stabilized spacecraft in a 250 km altitude circular orbit. The model is 
constructed as follows: 

i 

(a) Get a random angle (theta) in the X-Z plane measured with respect to the +Z 

direction (parallel to the spacecraft RAM-direction); 

• Get a new random number -1 < RAN#< 1 

• Compute the angle theta = invcos(RAN #) 

(b) Get a random angle (phi) in the X-Y plane measured with respect to the +X 

direction (perpendicular to the spacecraft RAM-WAKE axis); 

• Get a new random number 0<RAN#<1 

• Compute the angle phi = 2 n ■ RAN# 
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NOTE: this combination of (theta, phi) distributions ensures that the IDPs will arrive 
isotropically, i.e. there is an equal number of IDPs per unit solid angle arriving at the 
spacecraft from all directions. 

(c) Get a random IDP arrival velocity (Vp) by randomly sampling the Erickson IDP 
velocity distribution; 

• Get a new random number 0 < RAN # < 1 

• For RAN it < X m Solve F m (V P )~ RAN#=0 


NOTE: The IDP velocity distribution is constructed from 4 velocity distribution functions, 
f m {V P ) (m = 1..4), which fit the Erickson data in 4 velocity regimes, where; 

f-\{V P ) = Bq + a,V p + a^/p + a 3 Vp 
^2 iy p ) = exp(-b-|V P ) 

^3 i^p) — c 0 + C-\V p + c-^/ P + c 3 V P +C^Vp 

U{ v p) = d<i exp{~dy P ) 


F,(v P )=a 

F 2 (V„) = (-b 0 /b,)exf,{-by P ) 

F>(Vp) = c 0 V p+ ^ + ^ + Z& 
F 4 (Vp) = Ho/d 1 )exp(-d 1 l/p) 


c a V 5 p 


The constants aO, al, a2, a3, bO, bl, cO, cl, c2, c3, c4, dO, dl, vO, vl, v2, v3, v4 
necessary to compute the Erickson distribution are shown in Table 1, below. The 
functions F m (V P ) (m = 1..4) are the cumulative distribution functions, i.e. the finite 
integrals of the functions f m (V P ). The values X m (m = 1..4) are computed as follows; 



X 4 = 1 
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a 

b 

c 

d 

v [km/s] ' 

0 

-3 . 472e+00 

1 . 749e+00 

6 . 429e+00 

1.2374+01 

9.0 

1 

6 . 135e-01 

1 . 290e-01 

-4 . 364e-01 

1 . 290d-01 

19.5 

2 

-3 . 183e-02 

- 

1 . 104e-02 

- 

52.5 

3 

5 . 062e-04 

- 

-1 . 235e-04 

- 

64.5 

4 

- 

- 

5 . 144e-07 

- 

75.0 


Table 1. Constants for Erickson IDP Velocity Distribution Model 


The normalization constant, X c , is computed as follows: 

v, v, v 3 v 4 

x c = j f i( v ) dv+ J f 2 (v)dv+ jf 3 (v)dv+ j f 4 (v)dv 
V 0 Vi v 2 v 3 

(d) Compute the velocity ( V re j) relative to the spacecraft; for each arriving IDP: 

Vrei = a/(^sc + Vp cos(theta )) 2 + VpS\n 2 (theta) 

where Vsc = the spacecraft orbital velocity. 

(e) Compute the impact angle (gamma) relative to the spacecraft surface: 

(V S c + Vp cosftheta) 
gamma = /nvcosl — ^ ^ — * - 

V V ral 

(f) Determine the relative number of IDPs striking a specific surface: 

• Compute relative flux f = ^ rel cos (9 amma ) 

Vp 

• Get new random number 0 < RAN#< f ^ 

where f max = + ^ sc 

Vpmin 

• Remove all IDPs for RAN#>f 

NOTE that V p™, is the minimum arrival velocity (V P ) in the distribution selected from the 
Erickson distribution; in this case V Pmin = 1 1 .35 km/s. 
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Erickson 

Vrei-RAM 

V re t-WAKE 

gamma-RAM 

gamma-WAKE 

MIN 

11.4 

9.19 

5.46 

1.12 

5.88 

MAX 

74.6 

79.0 

70.7 

87.7 

88.9 

Points 

8192 

2512 

406 

2512 

406 

Mean 

20.2 

23.6 

18.3 

39.6 

53.1 

Median 

17.5 

21.7 

14.7 

38.7 

55.7 

RMS 

22.1 

25.0 

21.4 

43.7 

56.4 

Std Dev 

9.02 

8.02 

11.0 

18.5 

19.2 

Std Error 

0.100 

0.160 

0.548 

0.369 

0.955 

Skewness 

2.5 

2.5 

1.9 

0.26 

- 0.44 

Kurtosis 

8.3 

9.6 

4.2 

- 0.59 

- 0.63 


Table 2. Tabulation of statistical data for the isotropic IDP environment at 250 km altitude as 
described in the text. V re i = impact velocity; gamma = impact angle with respect to the specific 
surface normal. The relative fluxes for the RAM surface and the WAKE surface are F = 1.00 and 

F= 0.16, respectively. 

The results are shown in Figures 9-13. The main outputs of this modeling effort 
are the impact velocity - impact angle distributions for the RAM and WAKE surfaces 
(Figures 9 & 10), showing the accessible velocity-angle combinations for the two 
surfaces. Some combinations are not accessible due to the interaction of spacecraft 
orbital velocity and IDP arrival velocity. Effectively, the RAM surface aperture increases 
above 2 it steradians as a function of spacecraft velocity whilst the WAKE surface 
aperture decreases. In the case of the WAKE surface it is clear that if the IDP arrival 
velocity component parallel to the spacecraft velocity vector is less than the spacecraft 
velocity then the IDP cannot impact the WAKE surface. 

The impact velocity distributions are modified relative to the Erickson meteor 
arrival velocity distribution by the spacecraft velocity. The mean of the velocity 
distribution is shifted to higher velocities for the RAM-fadng surface, and shifted to lower 
velocities for the WAKE-facing surface (Figure 1 1 ). The dynamics of the system also 
flattens the WAKE distribution relative to the RAM distribution, spreading the impact 
velocities more widely. This is manifested by a reduction in the Kurtosis relative to the 
Erickson distribution. The statistics of the impact velocity and impact angle distributions 
are shown in Table 2. The Erickson meteor distribution is included for comparison. 

The Space Power Institute hypervelocity impact facility currently accelerates 
particles in the diameter range 50 - 400 pm to velocities in the range 5-13 km/s. Under 
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these conditions, the hypervelocity impact experiments reported here are representative 
of up to -30% of IDP impacts on WAKE-facing surfaces, but only -0.1% of RAM-fadng 
surface impacts (see Figure 13). For space debris particles, the average impact velocity 
in a circular orbit at 30° inclination is 9.8 km/s and therefore the facility can access up to 
-50% of the space debris impacts for RAM and WAKE surfaces combined. To sample 
these modified impact angle distributions impact experiments were conducted for a 
variety of impact angles varying from 30° to 75° and including a range of impact 
velocities from 5-13 km/s. 

4.2. Energy Dispersive X-ray Analysis 

To further understand the effects of impact obliquity and to provide support for the 
analysis of impact craters recovered from space-exposed quasi-infinite target surfaces, it 
is necessary to consider the abundance and distribution of impactor residues in the 
resultant craters. A program of impactor residue analysis using EDXS techniques was, 
therefore, executed. Samples on which hypervelocity craters were produced were 
placed in an SEM (Cambridge StereoScan 200) and the X-ray spectrum, derived from 
the excited states of the atoms present in the target/residue materials, was measured 
using an EG&G Ortec System 5000 X-ray spectrum analyzer. Compositional 
calculations were done using the SEM community-standard AutoZAP II program. 

For all analysis runs the accelerating voltage was set at 20 kV, which allowed the 
substrate (copper) signal to become significant. Elapsed live time was 100.00 seconds 
and the beam current varied from 0.3 nA to 4.6 nA. In all cases the specimen tilt angle 
was 45.0 degrees with respect to the X-axis and 0.0 degrees with respect to the Y-axis. 
Insertion depth was 60.0 mm with a typical working depth of between 1 5-1 9 mm. 

Quantitative data, corrected for substrate contamination (Si, Al), from which 
relative abundances of elements can be computed for each of the craters, are presented 
in Appendix A. Careful orientation of the crater with respect to the electron gun axis and 
the X-ray detector minimizes the geometric effects of deep craters. In section 4.2.1 a 
procedure for characterizing the effects of geometry is outlined to enable relative 
abundance computations to be made. 

4.2.1. Geometrical Effects 

The issue of electron beam and emission X-ray photon shadowing due to the geometry 
of the impact craters is quite significant. Essentially, the electron beam is incident on the 
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target surface which is inclined at 45° to the electron beam axis. The direction from the 
analysis site to the X-ray detector is perpendicular to the electron beam axis, and 
therefore under normal conditions where the target surface is flat the 45° target 
inclination maximizes the X-ray photon count. A problem arises when the target surface 
in the vicinity of the analysis site is not flat, i.e. in the case of an impact crater. By 
observing the X-ray spectra it becomes apparent that the relative Cu Ka and La line 
intensities vary with location dependent on the amount of electron beam and X-ray 
photon shadowing. Figures 14 through 18 show the X-ray spectra for the impact site 
E50-7-M2 (an 8 km/s impact at 45° incidence). It is clear from these spectra that the 
ratio of La to Ka line intensity varies significantly. Since the relative elemental 
abundance computations are dominated by the Ka counts and the Ka lines for Mg, Si, 
and Al are close in energy to the Cu La line it is necessary to re-scale the counts for Mg, 
Si, and Al to account for the shadowing effects. The rescaling was accomplished by 
determined the Cu La:Ka ratio for each of the analysis locations (1 - 5) within each 
crater. Figure 19 shows a plot of the average La:Ka peak height ratio as a function of 
location normalized to location (0) where location (0) is the original undamaged surface 
of the target. Location (1) includes location (1) and (5) data, i.e. on the crater lips. 
Location (2) includes location (2) and (4) data, i.e. on the crater interior wall. 

4.2.2. Substrate Contamination Effects 

The copper substrate is not 100% pure copper as it contains Si inclusions and Al 
contamination due to the polishing process. To correct for substrate contamination by Si 
and Al, EDXS analysis was conducted at undamaged target surface sites close to crater 
sites. Therefore, the relative composition of the substrate was determined for target E50 
to be: 


(Mg, Si, Fe, Al, Cu) = (0.00%, 4.41 %, 0.00%, 7.36%, 88.23%) 

i 

Corrections were made to separate the Si substrate and residue components to allow 
the computation of the residue:substrate ratio. An estimate of the Si substrate 
component was made for each analysis location by assuming that: 

Si substrate _ 4.41 
GUatomic% 88.23 

Then; 

Si residue = tom/c % — Si substrate 

where atomic% is the computed atomic percentage of the element at the analysis site. 
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4.2.3. Data 


Four preliminary experimental hypervelocity impact shots were executed, from which 
resulted in excess of 200 impact sites, identified using an optical scanning system. The 
minimum identifiable site diameter was -5 urn. The X-Y locations of all the sites were 
logged for subsequent correlation with streak record data. Shot D76 (normal incidence) 
was conducted to confirm gun performance and to determine the typical morphological 
structures that could be expected for olivine-copper impacts. 

Thirteen impact sites from the four shots were then selected for detailed analysis. 
These, listed below in Table 3, are representative of the numerous impact sites 
observed. Olivine residue was present in all of the selected craters, implying that 
residue retention can be expected at normal components of impact velocity up to 
12 km/s for the materials combination used in this set of experiments, and also for 
meteoroids striking aluminum targets if density ratio scaling is assumed to hold. 
Preliminary modeling of the flux dynamics of flat plates exposed to an isotropic 
meteoroid flux distribution (with the Erickson velocity distribution) indicates that at least 
12% of meteoroid craters on the RAM (East) surface should retain residue at levels 
detectable using EDXS, while at least 65% of WAKE (West) surface craters should 
retain similar residue. 


Site # 

Angle 

[deg] 

Velocity 

[km/s] 

Residue Appearance 

D76-01-16 

0 

4.8 

white, granular, crystalline 

D76-03-10 

0 

3.0 

white, granular, crystalline 

D76-03-13 

0 

9.9 

transparent, glassy 

D76-05-04 

0 

8.0 

white, granular, crystalline 

D76-05-05 

0 

9.6 

transparent, glassy 

D76-05-26 

0 

12.0 

transparent, glassy 

D76-05-32 

0 

12.0 

transparent, glassy 

D80-05-123 

45 

10.5 

transparent, glassy 

D78-05-23 

60 

8.5 

white, granular, crystalline 

D78-03-17 

60 

10.5 

transparent, glassy 

D79-01 -02 

75 

5.7 

white, granular, crystalline 

D79-03-04 

75 

12.5 

white, granular, crystalline 


Table 3. Preliminary sites selected for detailed SEM and EDXS analysis. 
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The normal incidence shot (D76) also provided a baseline for analyzing the 
residue structure within the craters. Two residue types manifested themselves in this set 
of craters: white, granular, crystalline; and transparent, smooth, glassy (amorphous). 
These morphologies are impact velocity dependent, the granular structures being 
evidence of lower shock intensity while the glassy residue results from shock-induced 
melting and re-solidification processes. 

A further experimental shot, E50, was executed and a detailed analysis of eleven 
(11) impact craters comprising five (5) EDXS analyses per impact crater was conducted. 
Rescaling of impact residue quantities was carried out in an attempt to account for the 
geometrical effects of the crater shape on the electron beam and emission X-ray photon 
count. The impact angle was set at 45° and the impact velocity spread was 3.7 km/s to 
12 km/s. The same two residue types manifested themselves as found in the 
preliminary impact crater set as shown in Table 4. 


Site it 

Angle 

[deg] 

Velocity 

[km/s] 

Residue Appearance 

E50-03-A 

45 

12.0 

transparent, glassy 

E50-06-B 

45 

11.0 

transparent, glassy 

E50-06-C 

45 

10.5 

transparent, glassy 

E50-07-AE 

45 

4.6 

white, granular, crystalline 

E50-08-I 

45 

8.5 

transparent, glassy 

E50-08-P 

45 

7.5 

white, granular, crystalline 

E50-08-N 

45 

7.9 

transparent, glassy 

E50-07-M1 

45 

8.0 

transparent, glassy 

E50-07-M2 

45 

8.0 

transparent, glassy 

E50-07-M3 

45 

8.0 

transparent, glassy 

E50-07-AB 

45 

3.7 

white, granular, crystalline 


Table 4. Sites selected for detailed SEM and EDXS analysis at 5 locations within each crater. 

The resultant data are given in Appendix A. In addition to the quantitative 
measurements taken at the five crater internal locations, a set of measurements were 
made at points on the undamaged surface of the target plates. These measurements of 
relative elemental abundances provide a baseline for the target material itself. 
Contamination of the copper substrate by both silicon and aluminum was detected and 
was corrected for in the computations of residue elemental abundances and ratio of 
residue to substrate at each analysis location. 
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5. 


DATA ANALYSIS 


5.1. Preliminary Impact Experiments 

Quantitative data, corrected for substrate contamination, for the preliminary impact 
experiments are shown in Table 5. Relative abundances of elements were then 
computed for each of the craters. Careful orientation of the crater with respect to the 
electron gun axis and the X-ray detector minimizes the geometric effects of deep craters. 
For normal incidence impacts location (1 ) is on the crater inner side wall and location (2) 
is on the crater lip. For oblique impacts (D78 and D80 sites) location (1) is situated on 
the leading edge inner side wall and location (2) is on the trailing edge inner side wall. 
For D79 sites, locations (2) and (3) are in downrange damage areas beyond the main 
crater trailing (downrange) edge. 


SITE# 

Mg 

Si 

Fe 

Cu 

Mg:Si 

RES% 

Yn 

mu 

D76-05-26(1 ) 

17.6 

5.4 

1.1 

75.9 

8.7 

21 

12.0 

n/a 

D76-01-17(1) 

28.8 

11.0 

2.0 

58.2 

3.4 

39 

2.9 

n/a 

D76-01-17(2) 

18.7 

5.4 

0.4 

75.5 

9.2 

22 

- 

- 

D78-03-17(1) 

23.7 

6.9 

2.4 

67.0 

6.1 

30 

5.3 

1.4±0.2 

D78-03-17(2) 

33.7 

10.9 

0.9 

54.5 

4.0 

43 

- 

- 

D78-05-23(1 ) 

25.0 

7.3 

0.9 

66.8 

5.8 

30 

4.3 

0.910.2 

D78-05-23(2) 

20.1 

8.2 

1.4 

70.3 

4.0 

27 

- 

- 

D79-03-04(1) 

30.1 

10.0 

1.7 

58.2 

4.1 

39 

3,2 

n/a 

D79-03-04(2) 

23.8 

8.1 

0.9 

67.2 

4.7 

30 

- 

n/a 

D79-03-04(3) 

18.5 

5.6 

0.5 

75.4 

8.3 

21 

- 

n/a 

D79-01-02(1) 

39.5 

17.1 

2.1 

41.3 

2.6 

57 

1,5 

n/a 

D79-01-02(2) 

23.2 

7.4 

0.9 

68.5 

5.3 

28 

- 

n/a 

D79-01 -02(3) 

20.4 

2.9 

0.5 

76.2 

n/a 

23 

- 

n/a 

D80-05- 123(1) 

15.9 

3.9 

1.3 

78.9 

41.8 

20 

7,4 

2.110.4 

D80-05-123(2) 

29.1 

10.4 

0.7 

59.8 

3.8 

41 

- 

- 


Table 5. Significant elemental composition by atomic percentage, Mg:Si ratio, and the residue 
percentage of the total sample. Substrate Si contamination counts were corrected for to obtain 
the Mg:Si ratio using the measured substrate atomic percentages i.e. Si=4.27%, Cu=95.73%. 

Two representative normal incidence impact sites were selected for more extensive 
morphological analysis. Site D76-01-17 was a 18413pm diameter crater, with a circular 
footprint and raised circumferential lips, exhibiting a p/D ratio of 0.42±0.02, computed 
from a measured depth of 77±3|xm. The impact velocity was 2.9±0.05 km/s. Coating the 
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interior of this crater was a whitish-green (when viewed optically) granular residue. 
Conversely, the lips of this crater were coated with a smooth, glassy residue, typical of 
that found in the interiors of the impact craters produced by faster impactors. 

At the other end of the velocity spectrum was site D76-05-26, produced by an 
olivine fragment traveling at 12.010.05 km/s. Again, the circular footprint and raised 
circumferential lips of a hypervelocity impact are evident. The p/D ratio was computed to 
be 0.4710.08 from a crater diameter (D) of 3812pm and a depth of 1812pm. The interior 
of this crater was coated with a smooth, glassy residue, transparent when viewed 
optically. 

The normal incidence crater set was characterized by a division into pre- 
dominantly granular residues and predominantly glassy residues occurring at an impact 
velocity of between 8.0 km/s and 9.6 km/s, which corresponds to a shock pressure of 
between 185 GPa and 240 GPa, using a computational scheme similar to that used by 
Ang (1990). The implications of this, for site D76-01-17, is that the impactor residue in 
the crater interior, in granular form, has undergone a significantly less intense shock 
metamorphosis than that deposited on the crater lips. 

The first trend to consider is that between the location (1) analyses, excluding 
D80-05-123 which results in significantly anomalous data. Normal component of impact 
velocity (V n ) correlates directly with residue Mg:Si ratio (Figure 20), implying shock- 
induced fractionation, i.e. the higher the velocity the more silicon is removed from the 
material matrix. This effect is confirmed by the appearance of the residue, i.e. the 
whitish-green granular residues are the least shock-fractionated (low Mg:Si ratio), 
whereas the smooth, glassy residues are the most shock-fractionated (high Mg:Si ratio). 

Also, residue percentage (RES%), determined by computing the ratio of residue 
(Mg, Si & Fe) counts to total counts (Mg, Si, Fe & Cu) for location (1) data shows 
correlation with normal component of velocity (V n ), with the highest velocities producing 
the lowest residue percentage, i.e. the residue layer thickness is apparently reduced for 
higher velocities (see Figure 21). A further consideration of the ratio of location (1) 
residue percentage to location (2) residue percentage for impact sites D78-03-17 and 
D80-05-123 indicates that for oblique angles of incidence more residue tends to be 
found on the trailing interior surface than the leading interior surface for craters with 
glassy, amorphous residues. The degree of downrange displacement appears to 
correlate with V n with higher normal components of impact velocity producing more 
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displacement. Consideration of the residue ratio (0.910.2) for site D78-05-23 shows that 
for sites with the white, granular, crystalline residue there is no apparent displacement of 
residue in the downrange direction. 

It is dear from these preliminary experiments that impact angle definitely affects 
the degree of shock-fractionation experienced by the impactor material and also the 
relative amounts of residue as a function of location. Effectively, two processes are at 
work: firstly, shock-fractionation; and, secondly, the displacement of the residue in the 
downrange direction as a function of the normal component of impact velocity. 

5.2. Detailed Analysis of Shot E50 

To provide further evidence of shock fractionation and residue distribution within impact 
craters, 1 1 craters (E50-series) were analyzed in detail. EDX spectra were obtained and 
the AutoZAP software used to compute the relative abundances of the major elements at 
each of five analysis locations within each crater. 


The Mg, Si, and Al atomic percentages were rescaled using the rescaling equation 
derived from the data shown in Figure 1 9 as follows: 


X'n - 


±n 

0.993-0. 267 W 


for X = Mg, Si, Al; A/ = 0..3 


where X' N is the rescaled atomic percentage. 

The Si and Al atomic percentages were then separated out into Si re sidue . S 'substrate, 
Al residue , and Al con tamination components using the procedure outlined in section 4.2.1 
above. The output results include: 


• the residue:substrate ratio as a function of normal component of impact velocity 
at the five analysis locations which allows the degree of residue retention within 
the crater to be determined; 

• the Mg:Si residue ratio as a function of normal component of impact velocity 
which allows the degree of shock fractionation of the residues to be determined 

• the ratio of downrange:uprange residues as a function of impact velocity which 
shows how the residue is distributed within the crater as a function of impact 
velocity. 
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Figures 22-26 show plots of residue:substrate ratio as a function of impact velocity 
normal component for shot E50 craters at analysis locations (1) through (5). All plots 
show a downward trend in retained residue as impact velocity increases. The 
downrange residue ratios, exemplified by Figures 22 and 23, show weak dependencies 
on impact velocity normal component whereas the uprange locations and the crater 
bottom show much stronger dependencies. 

Considering Figures 27-31 there appears to be correlation between shock fractionation 
and impact velocity normal component, although it exhibits a weak dependency. These 
data confirm the findings made during the analysis of the preliminary D-series shots. 
However, the spread in the data indicates that although shock fractionation is occurring 
to a greater or lesser extent, the quantification of such a process is far from being 
complete. It is also possible that shock fractionation effects are largely localized due to 
variations in shock compression release times. Such variations would result from the 
irregularity of the impactor in terms of shape, i.e. in these experiments and in space 
hypervelocity impactors are generally irregular in shape and so shock fractionation 
cannot be used as a tool for aiding in the determination of impact velocity. 

Figure 32, showing the plot of ratio of downrange residue percentage to uprange residue 
percentage as a function of impact velocity, is extremely interesting as it shows that as 
impact velocity increases the distribution of residues within the crater becomes more 
uniform. This is manifested by the convergence to an uprange:downrange residue ratio 
of unity. Such a convergence correlates with the tendency of impact crater geometry to 
converge to a hemispherically symmetric geometry as impact velocity increases. 

6. CONCLUSIONS 

The results from this work indicate that residue retention and residue distribution are 
definite functions of impact velocity and impact angle. It has been shown that residue 
percentages in the interior of impact craters depend heavily on the normal component of 
impact velocity. Likewise, there is also a tendency for residues to exhibit shock 
fractionation, again dependent upon the normal component of impact velocity. This 
latter phenomenon is not a pronounced as the residue retention process. 

The EDX analysis is accurate to ±0.01 percent under the conditions of the 
experimental analysis procedures used in this work. This implies that residues may still 
be detectable within impact craters at normal components of impact velocity up to 15 
km/s and perhaps even as high as 20 km/s. 
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One major finding has been the fact that as impact velocity increases the 
distribution of residues within the crater becomes more uniform as the 
uprange:downrange residue ratio converges to unity. This correlates with the 
hypervelocity impact phenomenon whereby craters become more hemispherically 
symmetrical as the impact velocity increases. 

The results obtained to date suggest that a combination of crater morphology and 
residue analysis can give evidence of impact azimuth even when the crater profile is 
near hemispherical as in the case of 30° - 45° impacts. The complex mechanics in 
oblique impacts where the residue is displaced downraoge and undergoes variations in 
shock fractionation as a function of position have begun to be unraveled. In particular, 
relative amounts and locations of residue do correlate with impact velocity and CSQ 
provide limits on that velocity. It is concluded that a combination of analysis techniques 
can yield further information on impact velocity, direction, and angle of incidence. It 
remains to be seen whether these findings can be used for LDEF crater analysis where 
the impact velocities are higher (15-25 km/s) and residues can only be detected by other 
techniques such as SIMS. 
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Figure 1 . Photomicrograph of an impact site formed by an olivine particle 
impacting with a velocity of 7.0 km/s at an angle of 30° with respect to the 
surface normal. Impact direction is indicated by the arrow. The interior of the 
crater is coated with olivine residue with visible particle grains. Note the 
circular footprint and symmetrical lips, [scale bar =100 pm] 


SPI 


NA G-1-1329 Final Report 


page 21 




Figure 2. Photomicrograph of impact site D80-05-1 23. Impact velocity was 10.5 km/s 
and angle of incidence was 45°. Impact direction is indicated by the arrow. Note the 
almost circular footprint and lack of asymmetry between leading edge and trailing edge 
lips. p/D ratio is 0.48, indicating an almost hemispherical crater which would be 
insdistinguishable from lower incidence impact craters. 


SPI 


NAG-1-1329 Final Report 


page 22 




particle residue 
[varying abundance] 


Figure 3. Schematic of idealized profile of normal-to-45° oblique incidence 
hypervelocity impact crater. The penetration depth (p) and crater radius (r) are related 
to the particle diameter (d) and normal component of velocity (v n ) similar to the 
relationship between penetration depth, particle diameter and impact velocity for normal 

incidence impacts. 



particle residue 
[varying abundance] 


Figure 4. Schematic of idealized profile of 45°-to-60° oblique hypervelocity impact 
crater. The minor ellipse (circle, in this case) parameters, penetration depth (p) and 
crater radius (r) are related to the particle diameter (d) and normal component of 
velocity (vn) similar to the relationship between penetration depth, particle diameter, and 
impact velocity for normal incidence impacts. The major ellipse parameters (a), (b) , 
and (Ap) are related to the angle of incidence and the shock propagation velocity in the 
target. 
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Figure 5. Schematic of the hypervelocity impact facility and the streak image velocity 

detection process. 
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Time t=0 is at the left-hand edge with time increasing to the right. The continous horizontal 
lines are the images of the non-pulsed LEDs, while the intermittent horizontal line is the image 

of the pulsed LED (T=300fis). 


particle direction of flight 



Figure 7. Typical crater top and side views showing locations of impact site quantitative 
EDX analyses. Location (1) is on the down range lip of the crater and location (5) is on 

the uprange lip. 
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number per BIN 


I DP arrival velocity relative to Earth 



Figure 8. Binned velocity distribution (Erickson meteor distribution) of IDPs entering 

near-Earth space. 
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RAM data 
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Figure 9. Plot of impact angle versus impact velocity for a RAM-facing surface of a 
gravity-gradient stabilized spacecraft in a circular orbit at 250 km altitude. The total flux 
number, N, is the flux over 4 k steradians. The flux, F, is the relative flux for the RAM- 

facing surface. 
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WAKE data 

N=8192; F=406; alt=250km 



impact velocity [km/s] 


Figure 10. Plot of impact angle versus impact velocity for a WAKE-facing surface of a 
gravity-gradient stabilized spacecraft in a circular orbit at 250 km altitude. The total flux 
number, N, is the flux over 4n steradians. The flux, F, is the relative flux for the WAKE- 

facing surface. 
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Figure 1 1. Velocity distributions for RAM- and WAKE-facing surfaces on a gravity- 
gradient stabilized spacecraft in a circular orbit at 250 km altitude. 


SPI 


NAG-1 -1329 Final Report 


page 30 



number per BIN 


N=8192; F = 2512; F =406 

RAM WAKE 

alt=250km 



BIN angle [deg] 


Figure 12. Impact angle distributions for RAM- and WAKE-facing surfaces on a gravity- 
gradient stabilized spacecraft in a circular orbit at 250 km altitude. 
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Figure 13. Probability distribution plot for impact velocity on both RAM and WAKE 
surfaces for a gravity-gradient stabilized spacecraft in a circular orbit at an altitude of 

250 km. 
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Figure 14. EDX spectrum for analysis location E50-7-M2-1. The La.Ka peak height 

ratio is 0.611. 
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Figure 15. EDX spectrum for analysis location E50-7-M2-2. The La.Ka peak height 

ratio is 0. 1 12. 


SPI 


NAG-1-1329 Final Report 


page 33 






Cfm 4K 




<-0.007 ksV Cursor - 0 C -0.007 ksV ) * O 19.196 ksV > 


Figure 16. EDX spectrum for analysis location E50-7-M2-3. The La:Kcc peak height 

ratio is 0.671. 


Cfs OK 
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Figure 1 7. EDX spectrum for analysis location E50-7-M2-4. The La:Ka peak height 

ratio is 0.557. 
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Figure 18. EDX spectrum for analysis location E50-7-M2-5. The La:Ka peak height 

ratio is 0.362. 



Figure 19. Normalized average La.Ka peak height ratio as a function of EDX analysis 
location. Location (0) is at the undamaged target surface where the electron beam and 
X-ray detector field of view is unobstructed by the crater geometry. Location (1) 
includes location (5) data, i.e. on the crater lip; location (2) includes location (4) data, 
i.e. on the crater interior wall. Location (3) is at the bottom of the crater 


SPI 


NAG-1-1329 Final Report 


page 35 








Figure 20. Graph of location (1) residues Mg:Si ratio showing the increase of shock 
fractionation with increasing normal component of impact velocity (V n ). This data is for 

the preliminary D-series shots. 



Figure 21. Graph of location (1) residue percentage versus impact velocity normal 
component (V n ). This data is for the preliminary D-series shots. 
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Figure 22. Plot of residue.substrate ratio as a function of impact velocity normal 
component for shot E50 craters at analysis location (1), i.e. on the downrange crater lip 

interior. 



Figure 23. Plot of residue.substrate ratio as a function of impact velocity normal 
component for shot E50 craters at analysis location (2), i.e. on the downrange crater 

interior wall. 
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Figure 24. Plot of residueisubstrate ratio as a function of impact velocity normal 
component for shot E50 craters at analysis location (3), i.e. in the bottom of the crater. 



impact velocity normal component [km/s] 


Figure 25. Plot of residueisubstrate ratio as a function of impact velocity normal 
component for shot E50 craters at analysis location (4), i.e. on the uprange crater 

interior wall. 
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Figure 26. Plot of residue.substrate ratio as a function of impact velocity normal 
component for shot E50 craters at analysis location (5), i.e. on the uprange crater lip 

interior . 



impact velocity normal component [km/s] 


Figure 27. Plot of residue Mg:Si ratio as a function of impact velocity normal 
component for shot E50 craters at analysis location (1), i.e. on the down range crater lip 

interior . 
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Figure 28. Plot of residue Mg:Si ratio as a function of impact velocity normal 
component for shot E50 craters at analysis location (2), i.e. on the down range crater 

interior wall. 



Figure 29. Plot of residue Mg.Si ratio as a function of impact velocity normal 
component for shot E50 craters at analysis location (3), i.e. in the bottom of the crater. 
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Figure 30. Plot of residue Mg:Si ratio as a function of impact velocity normal 
component for shot E50 craters at analysis location (4), i.e. on the uprange crater 

interior wall. 



impact velocity normal component [km/s] 


Figure 31. Plot of residue Mg:Si ratio as a function of impact velocity normal 
component for shot E50 craters at analysis location (5), i.e. on the uprange crater lip 

interior . 
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Figure 32 . Plot of ratio of downrange residue percentage to uprange residue 
percentage as a function of impact velocity. 
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APPENDIX A: Shot E50 EDXS Data 
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0.00 

4.57 


147.92 


0.459 


0.00 


53.88 


88.32 

86.92 

71,66 


0.192 


0.00 


41.09 


0.459 


115.29 


19.39 


8.89 


0.00 


0.726 


110.48 


20.51 


11.45 


6.29 


RES:SUB 


23.6 


13.1 


10.3 


0.55 


54.7 


36.4 


23.5 


1.41 


1.69 


39.2 


21.5 


27.6 


13.1 


0.61 


10.1 


2.73 


59.3 


16.8 


0.26 


75.4 


18.6 


10.4 


64.9 


0.35 


0.13 

0.00 

0.00 

0.00 

0.00 


-2.5 


. , 0.00 
- 0.00 
0.00 
0.00 
0.00 
























































































































































































































100.00 


100.00 


100.00 


100.00 


88.23 


13.09 


114.03 


48.43 


50.78 


12.19 


80.45 

69.54 


17.26 


SITE ID 


20.83 


E50-6-C1 


SITE ID 


115.81 


0.726 


44.59 


13.17 


129.80 


0.459 


45.10 


301.42 


0.192 


196.20 


53.07 


184.61 


0.459 


114.73 


41.68 


SITE ID 


E50-7-AE1 


RES:SUB 

Mg:Si 

0.92 

4.33 

0.68 

7.16 

4.63 

3.88 

5.74 

2.84 

1.43 

4.13 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


other 

total 

0.00 

108.51 

0.00 

187.89 

0.00 

0.00 

150.87 

121.25 

0.00 

111.22 


Fe 

Al 

Cu 

0.4 

0.0 

49.7 

0.8 

0.0 

56.7 

0.4 

0.0 

16.9 

1.1 

0.0 

14.1 

0.6 

0.0 

39.1 


SITE ID 

loc# 

SF 

E50-6-C1 

1 

0.726 

E50-6-C2 

2 

0.459 

E50-6-C3 

3 

0.192 

E50-6-C4 

4 

0.459 

E50-6C5 

5 

0.726 


















































































































































































































SUBSTRATE 


sub 1 






7.36 88.23 


totals 


100.00 


100.00 


100.00 


0.00 100.00 


SITE ID 




SF 


0.726 


0.459 


0.192 


0.459 


0.726 




Mq i Ma sea 


63.33 


45.34 


108.28 


99.43 


29.06 




Si Si sea 


20.41 


10.39 


27.03 


36.38 


16.67 


iy&iicif 



Al I Al sea 


0.00 


0.00 


0.00 


14.55 


11.46 




total 


122.95 


129.51 


209.34 


181.35 


100.37 


SITE ID 


E50-8-I1 

51.5 

16.6 

0.8 

0.0 

31.1 

1.6 

E50-8-I2 

35.0 

8.0 

1.4 j 

0.0 

55.6 

2.8 

E50-8-I3 

51.7 

12.9 

0.7 

o.o" 

34.6 

1.7 

E50-8-I4 

54.8 

20.1 

0.8 

8.0 

16.2 

0.8 

E50-8-I5 

29.0 

16.6 

1.1 

11.4"' 

” 42.0 

2.1 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


RES:SUB 


2.06 


0.71 


1.75 


4.07 


0.94 


SITE ID 


mmm 



SF 


0.726 


0.459 


0.192 


0.459 


0.726 



Ma I Ma sea 


57.47 


43.99 


156.46 


114.66 


46.02 



mm 


Si sea 


19.86 


10.17 


38.28 


42.81 


14.37 


Al Al 


llll * >11 





m 



total 


134.83 


127.94 


267.03 


201.87 


111.47 


SITE ID 

Mg 

Si 

fe 

Al 

Qj 

Si_sub 

Si_res 

E50-8-P1 

42.6 

14.7 

0.8 

5.4 

36.4 

1.8 

12.9 

E50-8-P2 

■Kina 

8.0 

1.0 

li II 

Kll 

K 

E50-8-P3 

58.6 

14.3 

0.5 

0.0 

26.6 

1.3 

13.0 

E50-8-P4 

56.8 

21.2 

0.8 

0.0 

21.2 

1.1 

20.1 

E50-8-P5 

41.3 

12.9 

0.9 

8.9 

36.0 

1.8 r 

11.1 


RES:SUB 


1.37 


2.58 | 

4.50 

3.50 ] 

2.82 

1.31 i 

3.72 
























































































































































































































SUBSTRATE 


sub 1 


0.00 4.41 


0.00 7.36 I 88.23 


SITE ID 


mmm 


SF 


0.726 


0.459 


0.192 


0.459 


0.726 


Ma sea 


46.64 


29.89 


89.58 


76.51 


50.33 


SITE ID 

■El 

Si 

Fe 

Al 

E50-8-N1 

39.2 

11.8 

0.6 

7.1 

E50-8-N2 

23.4 


3.1 

11.1 

E50-8-N3 

42.2 

10.2 

1.1 

13.0 

E50-8-N4 

46.0 

14.7 

0.6 

13.0 

E50-8-N5 

41.2 

13.0 

0.6 

12.4 


SITE ID 


mmm 




0.726 


0.459 


0.192 


0.459 


0.726 


mm 


Masca 


45.94 


60.15 


73.28 


51.46 


38.60 


SITE ID 


E50-7-M1-1 


m sam 


E50-7-M1-3 

39.2 

E50-7-M1 -4 

36.4 

E50-7-M1-5 

33.9 






















































































































































































































SUBSTRATE 


sub 1 


0.00 4.41 


0.00 7.36 88.23 


other totals 


100.00 


100.00 


100.00 


0.00 100.00 


SITE ID 


SF 


0.726 


0.459 


0.192 


0.459 


0.726 




46.93 


20.72 


168.33 


29.87 


17.60 


Si I Si sea 


14.53 


4.40 


53.85 


7.32 


5.32 


Al sea 


0.00 


0.00 


0.00 


0.00 


0.00 


Cu 1 other J total 


116.84 


113.59 


279.53 


120.11 


104.21 


SITE ID 

Kl 

Si 

Fe 

Al 

E50-7-M2-1 

40.2 

12.4 

0.7 

0.0 

E50-7-M2-2 

18.2 

3.9 

0.0 

0.0 

E50-7-M2-3 

60.2 

19.3 

0.4 

0.0 

E50-7-M2-4 

24.9 

6.1 

0.0 

0.0 

E50-7-M2-5 

16.9 

5.1 

0.0 

0.0 


Al con 


0.0 


0.0 


0.0 


0.0 


0.0 


RESrSUB 


1.04 3.98 


0.22 #### 


3.73 3.30 


0.38 9.41 


0.22 ! 13.66 


SITE ID 


SF 


0.726 


0.459 


0.192 


0.459 


0.726 


Ma sea 


23.25 


33.86 


80.94 


59.30 


14.26 


Si t Sisca 


6.71 


8.13 


18.70 


19.56 


2.95 


Al sea 


0.00 


0.00 


0.00 


10.85 


0.00 


total 


108.21 


122.71 


180.51 


148.44 


105.41 


SITE ID 


E50-7-M3-1 


\rnmm 


E50-7-M3-3 

44.8 

10.4 

0.0 

E50-7-M3-4 

40.0 

13.2 

0.6 

E50-7-M3-5 

13.5 

2.8 

6.6 



■■■ 

Al 

cu 

0.0 

72.3 



0.0 

44.8 

i warn 

39.0 

0.0 

83.6 


Sisub 


RES: SUB 


0.32 


0.4 




1.13 

5.52 

1.17 

3.56 

0.14 

-9.99 

























































































































































































































totals 


other 


Al 

Cu 

0.0 

28.2 

0.0 

13.5 


1.4 

18.1 

0.0 

0.0 

0.7 

25.8 

0.0 

0.0 

0.4 

27.3 

0.0 

0.0 

1.0 

19.7 

0.0 

0.0 

0.8 

25.4 

0.0 

0.0 



























































































































